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Abstract 

We report the first observation of CP violation in B° — > 7r + 7r~ decays based on 152 million 
T(4S) — > 55 decays collected with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. We reconstruct a 5° — > 7r + 7r~ CP eigenstate and identify the flavor of the accompanying 
B meson from its decay products. From the distribution of the time intervals between the two 
B meson decay points, we obtain A nn = +0.58 ± 0.15(stat) =L 0.07(syst) and S n7T = —1.00 ± 
0.21(stat) ± 0.07(syst). We rule out the CP-conserving case, = S n7T = 0, at a level of 5.2 
standard deviations. We also find evidence for direct CP violation with a significance at or greater 
than 3.2 standard deviations for any S n7r value. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the standard model (SM) of elementary particles, CP violation arises from the 
Kobayashi-Maskawa (KM) phase []J in the weak interaction quark-mixing matrix. In par- 
ticular, the SM predicts CP asymmetries in the time-dependent rates for B° and B° decays 
to a common CP eigenstate [2j]. Comparison between SM expectations and measurements 
in various CP eigenstates is important to test the KM model. The B° — > decay 0], 

which is dominated by the b — > uud transition, is of particular interest, and is sensitive to 
the CP- violating parameter 2 - Direct CP violation may also occur in this decay because 
of interference between the b — > u tree (T) and b — > d penguin (P) amplitudes pj]. 

In the decay chain T(45) — > B°B° — > (7r + 7r~)/ tag , where one of the B mesons decays at 
time to the CP eigenstate 7r + 7r~ and the other decays at time t ta , g to a final state /+ag 
that distinguishes between B° and B°, the decay rate has a time dependence given by |2j 

e -\At\/r B o 

V^(At) = — [1 + q ■ {5™ sin(Am d At) 

At b o 

+A 7T7T cos(Am d At)}} , (1) 

where t#o is the B° lifetime, Am^ is the mass difference between the two B° mass eigenstates, 
At = t n7T — t tag , and the 6-flavor charge q = +1 (—1) when the tagging B meson is a B° 
(B°). S nn and A nn are mixing-induced and direct CP- violating parameters, respectively. 

Belle's previous results for B° — > 7r + 7r~ [5], based on a 78 fb" 1 data sample (85 x 10 6 BB 
pairs), suggested large direct CP asymmetry and/or mixing- induced asymmetry while the 
result by the BaBar collaboration based on a sample of 88 x 10 6 BB pairs did not jfjj. In 
this Letter, we report a new measurement with an improved analysis that incorporates an 
additional 62 fb _1 for a total of 140 fb _1 (152 xlO 6 BB pairs) that confirms Belle's previous 
results with much greater significance. 

The data were collected with the Belle detector 0] at the KEKB asymmetric-energy e + e~ 
collider Q, which collides 8.0 GeV e~ and 3.5 GeV e + beams. The T(45) is produced with 
a Lorentz boost of f3j = 0.425 nearly along the electron beamline (z). Since the B° and 
B° mesons are approximately at rest in the T(4S) center-of-mass system (cms), At can 
be determined from Az, the displacement in z between the 7r + 7r~ and / tag decay vertices: 
At ~ {z n7T — z ta g) / file = Az/ftjc. The reconstruction method of the vertex positions remains 
unchanged from the previous publication (if . 

We use oppositely charged track pairs that are positively identified as pions to reconstruct 
B° — > 7r + 7r~ candidates. The pion efficiency is 91% and 10.4% of kaons are misidentified as 
pions. We select the B meson candidates using the energy difference AE = — E^ m 



and the beam-energy constrained mass M bc = y (E^ m ) 2 — (p ™) 2 , where E™* m is the cms 
beam energy, and E^ and p c ^ s are the cms energy and momentum of the B candidate. 
The signal region is defined as 5.271 GeV/c 2 < M bc < 5.287 GeV/c 2 and \AE\ < 0.064 
GeV, corresponding to ±3<r from the central values. To suppress the e + e~ — > qq continuum 
background (q = u, d, s, c), we form signal and background likelihood functions, Cs 
and Cbgi from the event topology, and impose requirements on the likelihood ratio LR = 
£s/(£s + £bg) f° r candidate events. We use the same event topology variables and the 
procedure that were used for the B(B° — > 7r°7r°) measurement [^|. 

The flavor of the accompanying B meson is identified from inclusive properties of particles 
that are not associated with the reconstructed B° — > 7r + 7r~ decay. We use two parameters, 
q [defined in Eq. (0)] and r, to represent the tagging information. The parameter r is an 
event-by-event, MC-determined flavor-tagging dilution factor that ranges from r = for no 
flavor discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort 
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data into six r intervals. The wrong tag fractions for the six r intervals, wi (I = 1,6), and 
differences between B° and B° decays, Awi, are determined from data [l0|. 

We optimize the expected sensitivity by using the improved likelihood ratio LR. We 
require LR > 0.86 for all r intervals. We include additional candidate events with lower 
signal likelihood ratio cuts (0.50, 0.45, 0.45, 0.45, 0.45, and 0.20) for different r intervals 
since the separation of continuum background from the B signal varies with r; we accept 
candidate events from twelve distinct regions in the LR-r plane. 




AE (GeV) 



FIG. 1: AE distribution in the M^ c signal region for B — > 7r + 7r candidates with LR > 0.86. 

Figure^shows the AE distribution for the B° — > 7r + 7r~ candidates that are in the M^ c sig- 
nal region with LR > 0.86 after flavor tagging and vertex reconstruction. In the Mt, c and AE 
signal region, we find 483 candidates with LR > 0.86 and 1046 candidates with LR < 0.86. 
The B° — ► 7r + 7r~ signal yield for LR > 0.86 is determined from an unbinned two-dimensional 
maximum likelihood fit to the M hc -AE distribution (5.20 GeV/c 2 < M bc < 5.30 GeV/c 2 
and —0.3 GeV < AE < 0.5 GeV ) with a Gaussian signal function plus contributions from 
misidentified B° — > K + tt~ events, three-body 5-decays, and continuum background. The 
fit yields 232±?g n + ix events and 82^ K + w events in the signal region, where the errors 
are statistical only. Extrapolating from the size of the continuum background in this fit, 
we expect 169 continuum events in the signal region. For LR < 0.86, the same procedure 
used in the previous publication [H[ yields 141 ± 12 n + ir~ events, 50 ± 8 K + tt~ events and 
855 continuum events in the signal region. The contribution from three-body S-decays is 
negligibly small in the signal region. 

The At resolution function R n7r for B° — ► 7t + tt~ signal events is formed by convolving four 
components: the detector resolutions for z n7T and z tag , the shift in the z tag vertex position 
due to secondary tracks originating from charmed particle decays, and the smearing due to 
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the kinematic approximation used to convert Az to At [10]. We assume — Rk-k and 
denote them collectively as R S i g . 

*A-7V7T *^7T7r 

are obtained from an unbinned maximum likelihood fit to the observed At 
distribution. The probability density function (PDF) for B° — > signal events (V% n ) is 

given by Eq.Q modified to incorporate the effect of incorrect flavor assignment. The PDF 
for B° — > K + ir~ background events is Vj^^At, wi, Awi) = ^-^-e~l At l/ r s° {1 — qAwi + q ■ (1 — 

2wi) ■ Ar-tt ■ cos(AmrfAt)}. We use Ak-k = as a default and include an effect of possible 
non-zero value for Ak-k in the systematic error. The PDF for continuum background events 
is V q q(At) = (1 + g-Ak g ){2^7e- |A * l/Tbk s + (1 - f T )5(At)}/2, where f T is the fraction of the 
background with effective lifetime r b k g and 8 is the Dirac delta function. We use ^4bkg — as 
a default. A fit to sideband events yields A^kg = 0.010 ±0.005. This uncertainty is included 
in the systematic error for A wn and S nn . All parameters of V qq (At) and R qq are determined 
from the events in the sideband region. 

We define the likelihood value for each (ith) event as a function of A n7T and S^: 

/+oo 
[{fZV^At^w^A^S^) 
-oo 

+f^VU^',w l )}-R sis (At l -At') 
+f™V qq {At') ■ R m {AU - At')]dAt' + f ol V ol {At t ). (2) 

Here the probability functions f™ (k = nn, Ktc or qq ) are determined on an event-by-event 
basis as functions of AE and M^ c for each LR-r interval (m= 1, 12) |5(. The small number of 
signal and background events that have large values of At are accommodated by the outlier 
PDF, V i, with fractional area / /. In the fit, S nn and A nn are the only free parameters 
determined by maximizing the likelihood function £ = Yl^Pi, where the product is over all 
B° — > 7r + 7r~ candidates. 

The unbinned maximum likelihood fit to the 1529 B° — > 7r + 7r~ candidates (801 B°- and 
728 i?°-tags), containing 372^11 n + ir~ signal events, yields A^ = +0.58 ± 0.15(stat) ± 
0.07(syst) and S W7r = —1.00 ± 0.21(stat) ± 0.07(syst). The correlation between A nn and 
S T7T is 0.286. As in our previous publication 0, we quote the rms values of the A nn and S. 
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distributions of the MC pseudo-experiments as the statistical errors of our measurement [11 
The usual fit errors from the likelihood functions, called the MINOS errors in the previous 
publication p, are In and for A n7T and S nn , respectively, in good agreement with 
the rms values above |l2| . In Figs. [2fa) and (b), we show the At distributions for the 264 
B°- and 219 -B°-tagged events in the subset of data with LR > 0.86. We define the raw 
asymmetry in each At bin by A = (N + — N-)/(N + + iV_), where iV+(_) is the number of 
observed candidates with q = +1(— 1). Figures Efc) and (d) show the raw asymmetries for 
two regions of the flavor-tagging parameter r. The effective tagging efficiency and signal 
purity is much larger in the 0.5 < r < 1.0 region. 

We test the goodness-of-fit from a x 2 comparison of the results of the unbinned fit and 
the At projections for B° -> tt+tt" candidates. We obtain xV D0F = 12.5/12 (7.6/12) for 
the At distribution of the B° (B°) tags. 

An ensemble of MC pseudo-experiments indicates a 26.7% probability of measuring CP 
violation at a level above the one we observe when the input values are A^ = +0.55 and 
S n7T = —0.84, which correspond to the values at the point of maximum likelihood in the 
physically allowed region (S% n + A\ v < 1); in this measurement it is located at the physical 
boundary (^ + 5^ = 1). 
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FIG. 2: The At distributions for the 483 B° -> vr+vr" candidates with LR > 0.86 in the signal 
region: (a) 264 candidates with q = +1, i.e. the tag side is identified as B°; (b) 219 candidates 
with q = — 1. (c) Asymmetry, A, in each At bin with < r < 0.5 and (d) with 0.5 < r < 1.0. The 
solid curves show the results of the unbinned maximum likelihood fit to the At distributions of the 
1529 B° -> 7T+7T- candidates. 

The systematic error is primarily due to uncertainties in the vertexing (±0.04 for Ann and 
±0.05 for S nn ) and the background fractions (±0.03 for A nn and ±0.02 for S W7r ). We include 
the effect of tag side interference [l3[ on ^4 7rw (±0.03) and 5 W7r (±0.01) in this analysis. Other 
sources of systematic error are uncertainties in the wrong tag fraction, physics parameters 
(Am d , t b o, and Ax n ), resolution function, background modeling, and fit bias. We add each 
contribution in quadrature to obtain the total systematic errors. The effect of the 3% charge 
asymmetry in the kaon misidentification rate is negligibly small. 

We perform a number of crosschecks. We measure the B° lifetime with the B — > 7r + 7r~ 
candidate events. The result, r^o = 1.46 ± 0.09 ps, is consistent with the world-average 
value A comparison of the event yields and At distributions for B°- and _B°-tagged 
events in the sideband region reveals no significant asymmetry. We select B° — > K + tt~ 
candidates by positively identifying the charged kaons. A fit to the 2358 candidates (1198 
signal events) yields Ak-k = —0.02 ± 0.08, in agreement with the counting analysis 
and Sk-k = 0.14 ± 0.11, which is consistent with zero. With the K + n~ event sample, we 
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determine t#o = 1.52 ± 0.06 ps and Am^ = 0.53lo;oy ps 1 , which are in agreement with the 
world average values |l4| . We check the measurement of A n7T using time- independent fits to 
the M hc -AE distributions for the B° and B° tags. We obtain A n7T = +0.73 ± 0.19, which is 
consistent with the time-dependent CP fit result. We also perform an independent analysis 
based on a binned maximum-likelihood fit to the At distribution. The result is consistent 
with that of the unbinned maximum-likelihood fit quoted here. 

The statistical significance of our measurement is determined from the same approach 
used in the previous publication p|. Figure El shows the resulting two-dimensional confidence 
regions in the A^ vs. «S 7r7r plane. The case that CP symmetry is conserved, A^ = S W7T = 0, 
is ruled out at the 99.999976% confidence level (CL), i.e., 1 — CL = 2.5 x 10~ 7 , equivalent 
to 5.2cr significance for Gaussian errors. The case of no direct CP violation, A n7T = 0, is 
also ruled out with a significance at or greater than 3.2<r for any S ww value. If the source of 
CP violation is only due to B-B mixing or AB = 2 transitions as in so-called superweak 
scenarios [l6j], then (S nn ,A n7T ) = (— sin20i,O). 1 — CL at this point is 8.4 x 10~ 4 , equivalent 
to 3.3a significance. 




FIG. 3: Confidence regions for A W7T and S^. The curves show the contours for 1-CL = 3.17 x 10 
(solid), 4.55 x 10~ 2 (dot-dashed), 2.70 x 10~ 3 (dotted), 6.34 x 10" 5 (dashed), and 5.96 x 10~ 7 (thick 
solid). 

Adopting the notation of Ref. the range of <p2 that corresponds to the 95.5% CL 
region for A^ and in Fig. El is 90° < <po< 146° for 0.15 < \P/T\ < 0.45 as used in the 
previous publication 5] and sin 20i = 0.736 1181. The result is in agreement with constraints 
on the unitarity triangle from other indirect measurements [HI. The 95.5% CL region for 
A^ and iS,,^ excludes \P/T\ < 0.17. 

In summary, we have performed a new measurement of CP violation parameters in 
B° -> 7T+7T- decays. We obtain A W7r = +0.58 ± 0.15(stat) ± 0.07(syst), and S W7r = 
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—1.00 ± 0.21(stat) ± 0.07(syst). We rule out the CP-conserving case, A nn = S nn = 0, at 
the 5.2cr level. We find evidence for direct CP violation with a significance at or greater 
than 3.2a. The constraints on 02 from our result are consistent with indirect measurements 
that assume the correctness of the SM. 
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